New constraints on Lorentz invariance violation from the neutron electric dipole 

moment 
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We propose an original test ol Lorentz invariance in the interaction between a particle spin and 
an electromagnetic field and report on a first measurement using ultracold neutrons. We used a 
high sensitivity neutron electric dipole moment (nEDM) spectrometer and searched for a direction 
dependence of a nEDM signal leading to a modulation of its magnitude at periods of 12 and 24 
hours. We constrain such a modulation to di2 < 15 x 10~ 25 ecm and ^24 < 10 x 10~ 25 ecm at 95 % 
C.L. The result translates into a limit on the energy scale for this type of Lorentz violation effect 
at the level of £ LV > 10 10 GeV. 



PACS numbers: 11.30.Cp, 11.30.Er, 13.40.Em, 12.60.-i 

The Standard Model of particle physics (SM) on the 
one hand and the theory of General Relativity on the 
other, are the two cornerstones on which our current un- 
derstanding of the Universe relies. Although of seemingly 
irreconcilable natures, the principle of Lorentz invariance 
is at the foundation of both theories. Unification of these 
two theories, including a consistent description of the four 
known interactions, is one of the main challenges of con- 
temporary physics. Among the many directions being 
explored, one of the most radical is to abandon space- 
time invariance under Lorentz transformations. A gen- 
eral framework to parameterize such Lorentz violating 
(LV) effects has recently been proposed It is based 
on the idea that diluted traces from primordial symmetry 
breaking can be observed via high precision experiments 
at low energies. 

Numerous such experiments have been performed over 
the last century. A first category of tests probes the 
photon sector, with a broad range of techniques from 
laboratory scale experiments to cosmological observa- 
tions 0. A second category deals with particles, includ- 



ing clock comparison experiments, spin polarized torsion 
pendula and accelerator based experiments. The cur- 
rent constraints on Lorentz violating vector and tensor 
background fields obtained from these experiments have 
recently been reviewed 0]. 

In this letter we report on an experimental limit for 
an interaction between a particle and an electromagnetic 
field resulting from a fundamental anisotropy of the uni- 
verse as recently proposed |4(. A nonrelativistic frame- 
work will first be developed followed by the description 
of the experimental procedure and the obtained results. 

Consider a nonrelativistic spin 1/2 particle in the pres- 
ence of electric and magnetic fields. Assuming rotational 
invariance, the form of the interaction potential is re- 
stricted to the simple form V = —fiaiBi — daiEi, when 
considering only the linear terms in the magnetic and 
electric fields Bi and J5j. Throughout this letter we 
adopt Einstein's repeated index convention and denote 
by <7j the Pauli matrices. Thus the interaction is de- 
scribed by only two quantities: the magnetic and electric 
dipole moments \i and d, respectively. Hence, allowing for 



2 



Lorentz violating background vector and tensor fields, in 
the spirit of [l[ and taking into account only linear terms 
in the electric and magnetic fields, the general form of 
the interaction potential becomes 

V = biCTi - (l,,rr,i:, - ^ija.Bj. (1) 

The first term bi is sometimes referred to as the cos- 
mic axial field, with the dimension of an energy. The 
most stringent limit is b < 10 -22 eV [5] but it has been 
searched for in numerous clock comparison experiments 
using different particles, including free neutrons @. The 
next terms and in Eq. (JlJ have the dimensions 
of an electric and magnetic dipole moment respectively. 
We will refer to dij (nij ) as the cosmic electric (magnetic) 
dipole tensor. They both violate rotation invariance be- 
cause they define privileged directions in the universe. 

The electric term leads to effects analogous to the 
electro-optical behavior of anisotropic media. If an elec- 
tric field is applied to a non centrosymmetric medium 
the latter becomes birefringent for light. This is known 
as the Pockels effect [3]. In our case, the vacuum itself 
is the medium and the particle spin corresponds to the 
polarization of the light. 

We probed these couplings by observing the spin 
precession of ultracold neutrons in the presence of a 
strong electric field and a weak magnetic field, using 
the RAL/Sussex/ILL spectrometer [8j dedicated to the 
search for the neutron electric dipole moment [9] . Under 
regular experimental conditions, a vertical Bq = 1 /iT 
magnetic field is applied parallel or antiparallel to a 
8.3 x 10 5 V/m electric field. The sensitivity to the elec- 
tric term is 10 4 times larger than to the magnetic one 
due to dimensional considerations. Thus, from now on 
we disregard the anisotropic magnetic moment and focus 
on the cosmic electric dipole tensor. 

The Ramsey's method of separated oscillating fields 
was used to measure the Larmor frequency of stored 
spin-polarized ultracold neutrons. Fluctuations of the 
magnetic field are corrected by means of a spin- 
polarized 199 Hg vapor as comagnetometer Both 
spin-polarized species (ultracold neutrons and mercury 
atoms) are stored in a cylindrical storage bottle (height 
h = 12 cm, radius r = 23.5 cm) during a measurement 
under vacuum conditions with a duration of 130 s. The 
storage bottle is composed of top and bottom electrodes 
coated with diamond-like carbon and of an insulating 
ring coated with deuterated polystyrene [ll| . The homo- 
geneous magnetic field Bq is generated by a coil inside a 
four-layer mu-metal magnetic shield. At the beginning of 
the precession time, transverse magnetic pulses are ap- 
plied to flip the polarization of both species by tt/2 onto 
a plane normal to Bq. The spin precession of mercury is 
monitored online by optical means. For neutrons a sec- 
ond coherent tt/2 pulse is applied at the end of the pre- 
cession time. The polarization is measured by sequential 



counting of the number of spin up and down neutrons 
leaving the storage volume. 

The cosmic electric dipole tensor has in general 9 com- 
ponents and it is convenient to split it into three parts: 

d^ = d°I ij +df j +d£, (2) 

where / is the identity matrix, d is the traceless, sym- 
metric part of the tensor and d A is the antisymmetric 
part. The first term d° is nothing else than the intrinsic 
EDM which actually does not violate rotational symme- 
try. The antisymmetric tensor is of rank 2 and dimension 
3, thus has 3 degrees of freedom. We define it as the axial 
vector 

d t = 2 € ijkdf kl (3) 

where e is the completely antisymmetric Levi-Civita ten- 
sor. Then the antisymmetric part of the interaction po- 
tential becomes 

V A = = (d A x E) • a. (4) 

This potential acts like a magnetic field orthogonal to the 
electric field. It is thus orthogonal to the main magnetic 
field Bq in the apparatus. To first order this additional 
field does not change the Larmor precession frequency 
and will therefore not be considered any further. This is 
not the case for the five terms arising from the symmetric 
part of the cosmic EDM tensor which are defined by 

(dxx dxY d X z\ 
d X Y dyy d Y Z , (5) 
dxz dyz dzz I 

with dzz = ~{dxx + dyy). These terms contribute to 
the Larmor precession frequency to first order. The Z 
axis is defined as the Earth rotation axis. While the 
Earth is rotating together with the vertical quantization 
axis and the applied electric field, these five contributions 
would show themselves as an EDM signal in three differ- 
ent ways: a steady shift d s tcady of the value of the EDM, 
a sidereal modulated part diA coming from the sidereal 
modulation of the direction of either the quantization 
axis or the electric field axis with respect to the static 
background tensor, and a part modulated at twice the 
frequency, di2, due to the combined effect of the mod- 
ulation of both axes. Taking into account the intrinsic 
EDM, we can write these three contributions as: 

dstcady = d Q + sin 2 Xd Z z + + dyy) (6) 



d 12 = cos 2 X^j -(dxx - dyy) 2 + d 2 XY (7) 

c?24 = 2 cos A sin \J d 2 xz + d YZ , (8) 
where A is the latitude of the experimental site. Then 

d(t) = d stea dy + dl2COs(2Q£ + ^i2) + d24COs(fii + 024) (9) 
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where SI = 27r/23.934 rad/hour is the sidereal frequency 
and (j>i2 and <p 2 4 are phases which depend on the defini- 
tion of the X and Y axes. 

Following the standard practice in the measurement of 
the neutron EDM using a comagnetometer, one considers 
the ratio R — / n //Hg ~ 30 Hz/8 Hz between the neutron 
and the mercury precession frequencies. In the presence 
of a homogeneous magnetic field and an electric field, this 
ratio depends on the direction of the latter according to: 



R(t) = 



7n 
7Hg 



2E 



d(t) 



(10) 



where 7„ and 7H g are the gyromagnetic ratios and d is the 
neutron EDM. We neglect a possible contribution from 
a time-dependent mercury EDM since the Hg nucleus is 
subject to Schiff screening of the electric field inside the 
atom 

We studied the time evolution of the correlation be- 
tween R and the electric field E during 5.6 days in 
December 2008 at the PF2 ultracold neutron beam- 
line at the Institut Laue-Langevin (ILL), Grenoble. An 
overview of the data is presented in Fig. [TJ where the 
variation AR of R around the mean value is plotted. 
While the main Bq field was pointing downwards, the 
electric field was reversed every 2 hours and some addi- 
tional data were taken without electric field to check for 
systematic effects. The statistical accuracy per cycle of 
the neutron frequency is given by [lOj : 
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2TrVNTa a e- T / T z 



30 ,uHz, 



(11) 



where typically N w 4600 is the number of neutrons 
per cycle, T = 130 s is the precession time and ctQ = 
0.86 ± 0.01 is the neutron polarization at the beginning 
of the precession time. The transverse neutron spin de- 
polarization time T 2 depends strongly on the magnetic 
field homogeneity, whereas the longitudinal depolariza- 
tion time Ti = 690 ± 80 s is attributed to depolarization 
occurring at wall collisions. The field homogeneity can 
be adjusted by a set of correction coils; the presented 
data was in fact taken in three different configurations 
for the currents in these coils. For the best magnetic 
field configuration, the value T 2 = 400 ± 38 s was ob- 
tained. In addition to the purely statistical error, we 
expect a fluctuation of the neutron frequency due to a 
random misalignment of the initial neutron spin after 
the mercury ir/2 pulse at the level of 20 /iHz in the worst 
case. Using only data at zero electric field, we indeed ob- 
serve a 18 /iHz non-statistical fluctuation. This error was 
added quadratically to the entire data set. The mercury 
cohabiting magnetometer was performing with a typical 
accuracy of 0.3 /iHz or 40 fT for an averaging time of 
130 s. Although negligible, the mercury contribution to 
the individual errors aR was taken into account. 

From a set of 1586 cycles, one can derive a value for 
^steady from the difference between R for the two different 
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FIG. 1: Variation of R as a function of time T for electric field 
up (upwards pointing red triangles) and down (downwards 
pointing blue triangles) and in the case of a null electric field 
(black dots). For each set of data, the mean value has been 
substracted. The data is folded modulo 24 h and then binned. 



directions of the electric field: 



Steady = (-3.4 ± 2.7 stat ) x 10"^ e cm. (12) 

Obviously, this <iteadv term is better constrained by 
the preceding work |9j using the same apparatus, i.e. 
d n < 2.9 x 10~ 26 ecm (90 % C.L.), where statistics has 
been accumulated for several years but where the time 
evolution has not been studied. However the fact that the 
present result was obtained in only about 5 days of data 
taking shows the high performance of the apparatus. 

Further, a Bayesian analysis was applied to the data 
to search for a time variation d(t), Eq. ©. First, the 
following Chi squared function is established: 
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X 2 (d 12 , 012,^24, 024 ) = ^2 



oR; 



(13) 
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where the sum runs over all data cycles and a 
Then the posterior probability density for di2 , d 24 is given 
by the likelihood function: 



L{d 12 ,d 24 ) = j-: 



cxp(-X 2 / 2 ) #12 



(14) 



where TV is a normalization coefficient. This function is 
plotted in Fig. [21 from which we deduce the following 
bounds: 

d 12 < 15xl0- 25 ecm (95% CL.) 



d 2i < 10 x 10- 



(95 % CL.) (15) 



This statistical limit could be affected by the main sys- 
tematic effect namely, a geometrical phase shift of the 
mercury precession frequency pro por tional to the elec- 
tric field and the vertical gradient [13|, |l4( : 
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FIG. 2: Isodensity lines for the posterior probability density 
function Eq. (|14[) . The probability inside the dashed (blue) 
line is 68 % and 95 % inside the solid (red) line. 



with ljq — |7H g i?o| the Larmor angular frequency and 
ojj. = 0.65 (vxy/r) 2 the effective radial velocity. In prin- 
ciple, modulations of the vertical gradient at periods of 
12 or 24 h would mimic the signal associated with new 
physics. The magnitude of the vertical gradients can be 
assessed from the value of Ro, the ratio of neutron to 
mercury precession frequency without electric field: 



Ro — 



fn 




-( 


/Hg 




7Hg V 



1 - 



dB-o/dz Ah 



(17) 



which originates from the vertical shift Ah of the neu- 
tron center of mass with respect to the mercury. Due 
to this gravitational effect, the neutrons and the mer- 
cury atoms do not average exactly the same magnetic 
field in the presence of a vertical gradient. By dedicated 
measurements of Rq and using Eqs. (JT6J) and (flTl) . it is 
possible to predict a false electric dipole moment signal: 
rffaisc = (1-2 ± 0.2) x 10~ 25 ecm. This shift in d stca dy is 
too small to be seen in the data with the given statis- 
tics. The uncertainty in df a i se has been calculated from 
the spread in the measured values Rq . These fluctuations 
are compatible with statistical fluctuations, in agreement 
with the previous measurement Q . In order to place an 
upper limit on the contribution of a modulated gradi- 
ent to our extracted limits di2 and c?24, one can take the 
uncertainty in Rq as the maximal amount of gradient 
fluctuations according to Eq. (fl"7|) . This then translates 
via Eq. (|16p into an upper limit of 2 x 10 -26 e cm as the 
systematic error in our limits due to gradient modula- 
tions. Given the current statistical sensitivity, this effect 
is negligible. 

Our result Eq. (|15l) can be simply interpreted on 
the basis of merely dimensional arguments. We de- 
note by £lv the energy scale associated with a viola- 
tion of Lorentz invariance. It is expected that d\i,di/± « 
ehcfEhv ■ This simple argument is supported by more so- 
phisticated arguments in a quantum field theory frame- 



work [J]. The limits in Eq. ([15]) correspond to a lower 
bound on the energy scale for Lorentz violation effects of 
10 10 GeV. This is far beyond energies accessible at parti- 
cle colliders (10 3 GeV) but still below the Grand Unifica- 
tion scale (10 16 GeV). Given that new physics is in gen- 
eral expected to be associated with a large energy scale, 
the proposed observables dy are indeed stringent tests of 
the Lorentz invariance complementary to the search for 
a cosmic axial field. 

A significantly improved sensitivity is expected in the 
near future with the same experimental installation, 
which has recently been moved to the Paul Scherrer Insti- 
tute. There it will benefit from a more intense ultracold 
neutron source [HI], and upgrades will allow for an even 
better control of systematic effects [1] . 
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